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ABSTRACT 


The  propagation  of  electromagnetic  radiation  through  the  atmosphere  is  a  crucial  aspect 
of  laser  target  acquisition  and  surveillance  systems  and  is  vital  to  the  effective  implementation 
of  some  Theater  Missile  Defense  systems.  Atmospheric  turbulence  degrades  the  image  or  laser 
beam  quality  along  an  optical  path.  During  the  past  decade,  the  U.S.  Air  Force’s  Geophysics 
Directorate  of  Phillips  Laboratory  collected  high  speed  differential  temperature  measurements 
of  the  atmt^heric  temperature  structure  parameter,  C^^  and  the  related  index  of  refiraction 
structure  parameter,  The  stratospheric  results  show  a  1-2  order  of  magnitude  increase  in 
day  turbulence  values  compared  to  night.  Resolving  whether  these  results  were  real  or  an 
artifact  of  solar  contamination  is  a  oitical  Theater  Missile  Defense  issue. 

This  thesis  analyzed  the  thermosonde  data  from  an  experimental  program  conducted  by 
the  Geophysics  Directorate  in  December  1990  and  found  strong  evidence  of  solar  induced 
artifacts  in  the  daytime  thermal  probe  data.  In  addition,  this  thesis  performed  a  theoretical 
analysis  of  the  thermal  response  versus  altitude  of  fine  wire  probes  being  used  in  a  new 
thermosonde  system  under  devefopment  at  the  Naval  Postgraduate  School.  Experimental  wind 
tunnel  measurements  were  conducted  to  validate  the  analytical  predictions. 
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I.  INTRODUCTION 


The  propagation  of  electromagnetic  radiation  through  the  atmosphere  is  a 
crucia)  aspect  of  many  laser  target  acquisition  and  surveillance  systems  and  is 
vital  to  the  effective  implementation  of  the  Theater  Missile  Defense  System  (TMD). 
Turbulence  degrades  the  quality  of  atmospheric  propagation  along  the  opticsd 
path  of  the  radiation.  Random  fluctuations  in  the  atmospheric  index  of  refraction 
produce  optical  turbulence  and  distort  the  phase  front  of  the  beam  [Ref.  1  ;p.  6-6}. 
The  index  of  refraction  structure  parameter,  C„*,  which  is  a  function  of 
atmospheric  wind  shear,  pressure,  and  temperature,  quantifies  optical  turbulence. 
Nearly  three  decades  of  measurements  show  that  C„‘  varies  both  horizontally, 
vertically,  and  with  time.  What  is  not  well-known  is  the  fine-scale  spatial 
distribution  and  the  intensity  of  turbulence  in  the  upper  troposphere  and  the 
stratosphere.  Atmospheric  temperature  fluctuations  characterized  by  the 
temperature  structure  parameter,  produce  index  cT  refraction  variations  by 
altering  the  air  density.  During  the  past  decade,  the  U.S.  Air  Force  Geophysics 
Directorate  of  Phillips  Laboratory  (PL/GP)  has  collected  measurements  of  this 
quantity  using  balloon-borne,  high-speed,  differential  temperature  thermosondes 
[Ref.  2:p.30].  Recently,  the  Naval  Postgraduate  School  (NPS)  has  used  both  fine 
thermistors  as  well  as  fine  wire  differential  thermocouples  to  measure  C/  in  both 
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day  and  night  conditions  [Ref.  3]  and  [Ref.  4],  The  GP  experiments  have  shown 
a  1-2  order  of  magnitude  increase  in  day  turbulence  values  compared  to  night. 
Attempts  to  explain  this  result  have  followed  two  lines  of  reasoning;  1)  either 
there  is  a  real  Increase  in  optical  atmospheric  turbulence  in  the  day,  or  2)  the 
data  (collected  by  two  3.8  pm  platinum-coated  tungsten  wire  probes  for  GP  and 
12.4  pm  chromel-constantan  thermocouples  for  the  NPS  Group)  are  unreliable 
in  the  day  due  to  systematic  effects  such  as  the  solar  heating  of  the  probes 
coupled  with  package  rotation,  pendulum  motion,  and  probe  vibration.  [Ref.  2;p. 
551 

This  thesis  investigated  two  systematic  error  effects  in  balloon  thermosonde 
systems  used  by  NPS  and  GP.  This  was  accomplished  in  two  phases.  The  first 
was  to  look  for  solar  induced  artifacts  in  the  daytime  thermal  probe  data.  The 
second  phase  was  to  predict  then  measure  the  thermal  probe  response  time  of 
the  NPS  probes  as  a  function  of  altitude.  TTie  initial  phase  of  this  investigation 
consisted  of  a  statistical  evaluation  of  a  December  1990  experiment  conducted 
by  GP  at  Holloman  Air  Force  Base,  New  Mexico,  that  focused  on  the  cross¬ 
correlation  of  measurements  taken  from  adjacent  pairs  of  probes  under  day  and 
night  conditions.  The  second  phase  was  an  analytical  development  and 
measurement  of  the  probe  response  time  versus  edtitude,  focusing  on 
thermodynamic  and  fluid  dynamic  factors.  These  calculations  were  compared 
with  wind  tunnel  measurements  to  validate  and  improve  current  probe  response 
models.  The  underlying  motivation  of  this  investigation  was  to  investigate 
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systematic  error  sources  in  critical  atmospheric  turbulence  measurements  needed 
for  effective  deployment  of  theater  missile  defense  systems. 
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II.  BACKGROUND 


A.  THEORY 

Lasers  emit  radiation  in  a  highly  directional,  collimated  beam  with  a  very  low 
angle  of  divergence.  The  I  earn  divergence  angle  B  can  be  "diffraction  limited," 
i.e.,  limited  by  the  wave  characteristics  of  light.  Under  vacuum  propagation 
conditions  the  divergence  angle  B  is 


where  ><  is  the  light  wavelength  and  0  is  the  diameter  of  the  aperture.  [Ref.  5:p. 
21]  Atmospheric  turbulence  degrades  tiie  spatial  coherence  of  a  beam  of 
electi’omagnetic  radiation.  The  Navier-Stokes  equations  govern  the  fine-scale 
random  fluctuations  of  turbulerit  flow  [Ref.  6:p.  262].  At  high  Reynolds  numbers, 
Pe~10®,  tfie  Navier-Stokes  equations  do  not  have  unique  solutions. 
Consequerttiy  we  must  use  statistical  models  based  on  semi-empirical, 
dimensional  analysis  and  physical  reasoning.  These  focus  on  the  mean  flow 
properties  of  frie  fluid  and  on  second  order  statistics  of  tite  fluctuations 
(correlations,  spectra,  and  structure  functions)  [Ref.  7;  p.  311].  Aooilimatedlaser 
beam  irtitiaily  propagates  with  a  uniform  wave  front  until  it  encounters  a  tuitjuierrt 
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medium  (the  atmosphere).  Atmospheric  temperature  inhomogeneities  cause 
random  fluctuations  in  the  atmospheric  index  of  refraction,  n,  [Ref.  1  :p.  6-6]  which 
distort  the  phase  of  the  wave.  The  result  of  this  distortion  on  an  initially 
collimated  laser  beam  is  to  induce  centroid  wander,  broadening,  and  break-up 
(in  strong  turbulence)  which  are  caused  by  a  general  reduction  in  the  spatial 
coherence  of  the  original  beam  (Rgure  1). 
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B.  STRUCTURE  FUNCTIONS  AND  PARAMETERS 


The  mean  square  difference  between  two  arbitrary  measurements  X  such 
as  temperature  is  defined  as  a  structure  function 

where  0,^  depends  on  the  separation  distance  [Ref.  8:p.  19].  If  we  assume 
homogeneous,  isotropic  turbulence,  then  Kolmogorov  and  Tatarski  [Ref.  8:p.  46] 
have  shown  that 

(3) 


where 


r»|ri>r2i=ri2. 


(4) 


This  dependence  depicts  Kolmogorov  turbulence  represented  by  the  structure 
parameter  C/. 

The  Index  of  refraction  stnjcture  parameter,  C„®  describes  the  fluctuations 
in  the  index  of  refraction  over  a  given  distance.  It  is  the  mean-square  statistical 
average  of  the  difference  in  the  index  of  refraction,  n,  between  two  points 
separated  by  a  distance  r^j  [Ref.  9;p.  397] 


■12 


2/3 


(5) 
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Similarly,  the  temperature  structure  parameter,  C/  is 


c 

'-T 


■12 


2/1 


(C) 


where  and  T,  are  the  ambient  atmospheric  temperatures  measured  along  a 
horizontal  path  of  length  r.  According  to  Kolmogorov,  this  relation  is  only  valid 
in  the  inertial  sub  range  described  by  lo<<r<  <L,,  where  1^  and  depend  on  the 
kinetic  energy  and  boundary  conditions.  The  inner  scale,  l^,  where  energy 
dissipation  by  viscosity  and  diffusion  takes  place,  is  typically  on  the  order  of 
millimeters  near  the  ground  and  perhaps  centimeters  in  the  troposphere.  The 
outer  scale,  L^,  the  largest  scale  of  isotropy,  is  usually  on  the  order  of  meters  or 
10s  of  meters  IRef.  1:p.  6-9]  and  depends  on  ttie  boundary  conditions. 

The  index  of  refraction,  n,  is  a  function  of  atmospheric  pressure,  humidity, 
temperature,  and  c^tical  wavelength.  We  may  neglect  the  effects  of  humidity  in 
all  but  the  most  humid  atmospheric  environments  and  can  disregard  the 
wavelength  dependence  at  optical  wavelengths.  For  optical/IR  wavelengths,  a 
simplified  form  of  this  critical  parameter  is 

n-l»79xio**-|,  (7) 
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where  P  is  atmospheric  pressure  (millibars)  and  T  is  the  atmospheric  temperature 
(Kelvin)  [Ref.  10:p,  232].  Taking  the  partial  derivative  of  Eq.  (7)  and  assuming 
isobaric  turbulence,  depends  on  Ct*  by 

C„2*(79x10'‘-^)2c^2.  (S) 

The  effects  that  turbulence  produces  on  a  laser  beam  are  seen  by 
measuring  the  amplitude  or  phase  fluctuations  of  the  final  image  at  ^e  target. 
A  common  measure  of  image  quality  is  the  Strehl  ratio  which  is  the  relative 
intensity  of  the  point  spread  function  measured  at  the  mean  centroid  location 
relative  to  an  ideal,  diffraction  limited  optical  system.  The  Strehl  ratio  depends 
on  the  spatial  coherence  of  the  laser  radiation  after  passing  through  the 
atmosphere  which  in  turn  depends  on  the  atmospheric  modulation  transfer 
function.  The  atmospheric  modulation  transfer  function  (MTF)  characterized  by 
the  transverse  coherent  path  length,  rg,  is  a  measure  of  the  electric  field 
transverse  autocorrelation  lengtfi,  the  distance  over  which  the  electromagnetic 
field  remains  coherent  (Ref.  1;p.  6-29].  Fried  showed  that,  for  plane  waves, 

ro=2.1(1.46lc2^C„M^)d^]'’/^  (9) 

0 

where  k  is  ^e  opticsd  wave  number  (2/rM)  and  L  is  the  optical  path  length  [Ref. 
9;p.  397].  Typical  values  of  r^  range  from  centimeters  in  regions  of  high 
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turbulence,  to  tens  of  centimeters  at  good  astronomical  sites,  to  meters  and 
perhaps  larger  in  low  turbulence  regions. 

C.  CURRENT  APPUCATIONS 

Some  alternatives  to  the  TMD  scenario  involve  near  horizontal  atmospheric 
propagation  at  altitudes  of  10-20  km.  One  initiative,  the  Airborne  Laser  Program 
(ABL),  involves  mounting  a  laser  in  an  aircraft  flown  above  the  tropopause  into 
the  stratosphere  to  avoid  clouds,  sub-visual  cirrus  (ice  crystals),  and  tropospheric 
turbulence  [Ref.  11].  The  upper  altitude  limit  arises  from  aircraft  attitude 
constraints  although  other  platforms  such  as  low  orbit  satellites  are  also  being 
considered  (space^sased  alternatives  are  subject  to  more  restrictive  political 
constraints).  The  lower  altitude  limit  essentially  is  the  bottom  of  the  stratosphere. 
The  reason  for  this  is  that  beam  propagation  in  the  troposphere  encounters 
higher  turbulence,  particularly  near  the  jet  sti'eam,  and  has  higher  absorption  and 
scattering  from  atmospheric  aerosols,  water  vapor,  ice  crystals,  and  clouds. 

With  adequate  phase  information,  adaptive  optics  can  pre-distort  the 
wavefront  of  tiie  beam  to  correct  for  turbulence  that  is  along  the  propagation 
path.  This  is  necessary  to  achieve  a  30  cm  -  1.0  m  diameter  spot  size  on  a 
target  at  a  range  of  400  km  (0.8  Arrad).  The  Airborne  Laser  Program  needs  an 
accurate  assessment  of  the  atmospheric  turbulence  along  stratospheric  optical 
paths  to  complete  a  systems  analysis.  This  was  the  driving  force  for  this 
investigation. 
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D.  PREVIOUS  RESULTS 


GP  and  NPS  have  conducted  numerous  experiments  to  measure  the 
temperature  structure  parameter  using  bailoon-bome  thermosondes  ^d  single¬ 
point  temperature  probes.  Thermosondes  measure  differential  temperature 
fluctuations  across  a  horizontal  path,  commonly  one  meter  [Ref.  12:p.  14].  The 
mean  square  fluctuation  in  temperature  between  two  3.8  //m  platinum-coated, 
tungsten  probes  separated  by  one  meter  gives  a  measurement  of  Cj  which,  in 
ttjrn,  gives  using  Eq.  (8)  [Ref.  13:p.  2].  The  NPS  single-point  temperature 
probes  measure  temperature  using  100  pm  thermistors  along  both  the  ascending 
and  descending  trajectories  of  the  balloon. 

Previous  experiments  and  analyses  have  found  a  clear  increase  in  the  index 
of  refraction  structure  parameter  above  -*8  km  during  daylight  hours.  In 
June  1990,  GP  conducted  a  series  of  thermosonde  launches,  in  conjunction  with 
NPS  optical  measurements.  Rgures  2  and  3  show  plots  of  from  both  day 
and  night  launches.  The  raw  data  minima  in  Figure  2,  a  night  flight,  approach 
the  sonde  noise  level  (solid  curve)  for  much  of  the  flight.  The  raw  data  in  the 
daytime  flight  (Figure  3),  almost  never  approach  the  sonde  noise  level  [Ref.  3:p. 
55].  This  anomaly  has  been  characteristic  of  all  day  flights  and  is  suspicious. 
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ALTrniDS  (km) 


Figurs  3.  Log  (Raw  C,*)  Profile  Collected  by  the  Geophys ics 
Directorate  on  Maul,  3  June  1990  (Day)  [Ref.  3:p.55]. 
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GP  and  Sparta,  Inc.,  published  a  scientific  report  in  August  1990  addressing 
temperature  perturbations  from  solar  heating  of  thin  layers  of  aerosol  and  ozone. 
They  found  that  there  was  some  evidence  to  support  the  postulate  of  temperature 
perturbations  due  to  solar  heating  of  ozone,  although  this  could  not  be  proven 
due  to  the  low  vertical  resolution  nature  of  O,.  GP  concluded  that  it  was  not 
probable  that  variations  in  C/  were  due  to  perturbations  in  background 
stratospheric  aerosols.  (Ref.  13:p.  14]  GP  also  investigated  the  possibility  that 
solar  radiation  heating  of  the  probe  support  wires  caused  instrumental  effects  but 
concluded  that  this  was  not  responsible  for  the  prevalent  diurnal  variation  [Ref. 
2:p.  55]. 

Previous  work  at  NPS  suggested  that  solar  heating  and  probe  wire 
vibrations,  caused  by  the  wind  from  the  balloon’s  ascent  coupled  with  the  probe 
package’s  oscillating  motion,  produced  an  AC  signal  that  seriously  degraded  the 
daytime  data  [Ref.  3].  An  actual  increase  in  daytime  turbulence  was  extremely 
unlikely  since  the  day-night  change  In  turbulence  occurred  within  minutes  of 
sunrise  or  sunset,  it  is  unlikely  that  this  rapid  onset  and  decay  is  a  real 
phenomena  because  large-scale  molecular  diffusion  of  energy  throughout  the 
abnosphere  should  not  be  possible  on  such  small  time  scales.  Optical 
measurements  collected  in  July  of  1990  showed  no  indication  of  a  change  in 
isoplanatic  angles  between  day  and  night  v^ereas  thermosonde  results  collected 
during  the  same  period  showed  persistent  day-night  differences.  Optical 
isopianabc  angles  are  sensitive  to  weighted  by  altitude  z’^  [Ref.  3:p.  9], 
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Consequently,  turbulence  in  the  5-15  km  region  has  the  dominant  contribution. 
Weitekamp  showed  that  simulated  solar  heating  of  a  GP  3.8  tungsten  probe 
In  the  laboratory  produced  a  0.1  to  0.2  C  temperature  increase.  This  was 
sufficient  to  account  for  the  obsen^ed  increase  in  the  day  values  when 
coupled  with  probe  vibration  and  sonde  rotary  and  pendulum  motion. 
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ill.  DIURNAL  VARIATION 


A.  DIURNAL  VARIATION  EXPERIMENT  2 

1.  Purpose 

The  Geophysics  Directorate  group  conducted  a  series  of  thermosonde 
launches  from  4  to  9  December  1990  at  Holloman  Air  Force  Base,  New  Mexico, 
to  study  the  diurnal  variation  in  turbulence  measurements  of  C„*.  The  purpose 
was  to  validate  or  disprove  the  existence  of  measurement  artifacts  that  were 
responsible  for  the  diumai  variations  observed  by  both  NPS  and  GP  in  previous 
experiments.  [Ref.  14]  These  artifacts  include: 


1.  Balloon  thermal  wake  in  combination  with  payload  motion  (pendulum 
and  rotation). 

2.  Differential  solar  heating  of  the  probes  due  to  misalignment  in 
combination  with  payload  motion  (pendulum  and  rotation). 

3.  Differential  solar  heating  of  the  probes  due  to  different  wire  geometry  or 
emissivity/reflectivity  in  combination  with  payload  motion  (pendulum  and 
rotation). 

4.  Differential  solar  heating  of  the  probes  due  to  probe  wire  oscillation, 
possibly  exacerbated  by  payload  motion  (pendulum  and  rotation). 


2.  Experimental  Procedure 

The  GP  thermosonde  system  measured  the  temperature  difference 
between  two  3,8  /jm  platinum-coated  tungsten  wire  probes  exposed  to  the 
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ambient  atmosphere.  The  probe  wires  were  flexible  s-shaped  or  u-shaped  wires 
approximately  0.45  cm  long  connected  to  two  probe  supports.  The  probes 
formed  elements  of  an  AC  Wheatstone  bridge  which  measured  the  difference  in 
resistance  between  them.  High  and  low  pass  filters  restricted  the  system 
bandwidth  to  the  range  .2  Hz  to  1000  Hz  at  the  -3dB  level.  The  thermosonde 
passed  the  temperature  difference  Tj  -  T,  through  a  root  mean  square  (RMS) 
module  with  a  4-second  time  constant  to  ground-based  monitoring  equipment 
every  four  seconds  which  logged  the  data  and  calculated  [Ref.  12;p.  15]. 
The  sondes  ascended  at  an  average  rate  of  approximately  5  meters  per  second. 
An  ascent  to  30  km  required  approximately  95  minutes.  Historically,  a  GP 
thermosonde  launch  consisted  of  a  thenmosonde  package  with  a  single  pair  of 
probes  with  a  1  m  spacing  connected  to  a  VIZ  MICROSONDE  1  telemetry 
package.  The  December  1990  launches  were  unique  since  they  used  two 
separate  thermosonde  packages  strapped  together  with  a  single  VIZ  telemetry 
package.  The  two  pairs  of  probes  were  arranged  in  1  m  differential  pairs  on  a 
single  50  x  50  mm  styrofoam  boom,  so  that  each  sensor  element  was 
approwmately  8  cm  from  the  corresponding  member  of  the  other  pair.  There 
were  two  probe  configurations  for  the  experiment.  On  most  flights,  both  peurs  of 
probes  were  parallel.  One  flight  had  one  pair  of  probes  parallel  and  the  second 
set  perpendicular  to  each  other. 


3.  Experimental  Conditions 

In  all,  five  successful  balloon  launches  occurred  during  December 
1990.  Four  were  day  launches  and  one  was  a  night  launch.  Local  sunset 
occurred  around  1653  hours.  Difficulties  in  controlling  the  heavier  dropsonde 
packages  precluded  any  data  taking  during  descent.  Table  I  contains  a 
summary  of  the  individual  sonde  launches. 

a.  4  December  1990 

Time  of  launch  of  DVE9465  was  1450  hours  local  and  both  sets 
of  probes  were  parallel.  The  last  good  data  were  obtained  at  29.03  km. 

b.  5  December  t990 

The  time  of  launch  of  DVE9460  was  1407  hours  local  on  5 
December  1990  and  both  sets  of  probes  were  parallel.  Good  data  were  obtained 
up  to  26.48  km. 

c.  6  December  1990 

The  time  of  launch  of  DVE9454  was  1330  hours  local  on  6 
December  1990.  Both  sets  of  probes  were  oriented  in  a  parallel  fashion.  The 
maximum  altitude  providing  useful  data  was  20.5  km. 

d.  7  December  t990 

The  time  of  launch  of  DVE9448  was  1357  hours  local  on  7 
December  1990.  One  set  of  probes  was  in  a  parallel  configuration  while  the 
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second  set  was  peipendicular.  This  was  the  only  launch  in  \AHiich  this  was  the 
case.  The  maximum  altitude  that  provided  useful  data  was  28.86  km. 

e.  8  Decemlier  1890 

DVE9462  was  the  only  night  launch.  Time  of  launch  was  1837 
hours  local.  Both  sets  of  probes  were  parallel.  Good  data  were  taken  up  to 
30.02  km. 


TABLE  I.  Diurnal  Variation  Experimervt  2  Launch  Summary,  Conducted  by  the 
Geophysics  Directorate  at  Holloman,  Air  Force  Base,  NM  from  4-8  December 
1990. 


DATE 

4DEC90 

5DEC  90 

6DEC  90 

7DEC  90 

8DEC90 

FUGHT 

DVE9465 

DVF9460 

DVE9454 

DVE9446 

DVE9462 

TIME 

(Local) 

1450 

1407 

1330 

1357 

1837 

PROBES 

par/par 

par/par 

par/par 

par/perp 

par/par 

MAX  ALT 
(km) 

29.03 

26.48 

20.5 

28.86 

30.02 
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B.  RESULTS  AND  ANALYSIS  OF  DIURNAL  VARIATION  EXPERIMENT  2 


1.  Data  Raduction 

The  telemetry  package  continuously  transmitted  data,  at  four  second 
intervals,  to  the  ground-based  monitoring  equipment,  resulting  in  data  sets  of 
approximately  1450  lines  for  a  90  minute  flight.  The  data  were  a  succession  of 
independent  4  second,  RMS  temperature  differences  in  the  form 

[<(r2-ri)2>]i/2  (10) 

Eq.  (6),  the  temperature  structure  function  equation,  along  with  the  temperature 
difference  data  and  the  r=1  m  spacing  of  the  probes  gave  a  C/  value  for  each 
altitude  at  which  measurements  were  made,  approdmately  every  20  m. 
Combining  these  data  with  measured  pressure  and  temperature  yielded 
values  at  each  altitude. 

2.  Statistical  Methodology 

To  determine  if  any  measurement  artifacts  contaminated  the  data,  we 
performed  three  statistical  analyses  of  the  data  taken  from  the  two  pairs  of 
probes.  The  first  involved  examining  plots  of  the  individual  data  points  vs. 
altitude.  This  can  reveal  major  inconsistencies.  The  second  involved  an 
examination  of  x-y  cross-correlation  plots  between  data  from  pairs  of  probes. 
The  third  involved  computation  of  the  correlation  coefficient  between  probe  pairs. 
The  temperature  differences  measured  at  each  successive  altitude  should  have 


a  high  correlation,  since  both  sets  of  probes  were  essentially  passing  through  the 
same  air  mass.  If  there  were  independent  measurement  artifacts,  however; 
inconsistencies  in  the  x-y  cross-correlation  plots  and  the  correlation  between  data 
could  possibly  provide  some  insight  as  to  the  source  of  the  error. 

This  thesis  used  the  Microsoft  EXCEL  spreadsheet  to  compute  the 
statistical  correlation  between  both  pairs  of  RMS  temperature  fluctuations.  The 
correlation  coefficient,  r,  was  computed  according  to  the  standard  statistical 
expression 

r  ( CT, .  CTJ  =  .  (11) 

V®CM®CT2 


where 


(Cy-rT)  , 

*X  **1  ^2 


(12) 


and  CT1  was  the  RMS  temperature  difference  for  one  pair  of  probes  and  CT2 
was  the  difference  for  the  second  pair.  Values  of  r^- 1  indicate  a  high  degree  of 
correlation;  values  of  r~0  indicate  little  or  no  correlation. 

For  analyses,  the  30  km  data  sets  were  divided  into  5  km  segments. 
We  also  computed  a  single  r  coefficient  for  each  flight  that  contained  edl  the  data 
points.  Applying  Eq.  (12)  to  each  vertical  5  km  region  (or  to  the  complete  flight 
in  the  case  of  the  total  data  set)  gave  the  standard  deviation  needed  in  Eq.  (11), 
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which  in  turn,  yielded  the  correlation  coefficient  between  each  pair  of 
measurements. 

The  measured  noise  level  of  the  instrument  packages  was  .002K 
Consequently,  (AT),^  is  undefined  for  values  less  than  .002K  [Ref.  12:p.  14],  All 
data  taken  that  produced  a  reading  of  :s.CX)2K  were  discarded  as  being 
statistically  invalid  because  it  would  not  be  possible  to  discriminate  between  the 
signal  and  the  noise. 

3.  Results 

a.  Diurnal  Variation 

Plots  of  Log  (Raw  vs.  altitude  show  the  typical  order  of 
magnitude  increase  in  optical  turbulence  in  the  day  compared  to  night  in  tiie  1 0* 
20  km  region.  Figures  4  and  5  of  the  7  December  flight  represent  a  daytime 
flight  and  are  also  the  only  plots  obtained  using  the  parailei/perpendicular 
configuration  (Appendix  A  contains  the  otiier  daytime  plots).  Rgures  6  and  7 
show  the  data  taken  by  the  two  sets  of  probes  on  the  8  December  night  flight. 


Figux*  4.  Log  (Raw  C.^)  vsT  Altitude  Profile  ^ytime) 
collected  by  the  Geophysics  Directorate,  from  1.2  to  30  km 
Holloman  Air  Force  Base,  NM,  7  December  1990  (Probe  Set  1) . 
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Figur*  5.  Log  (Raw  C/)  vs^  Altitude  Profile  (Daytime) , 
Collected  by  the  Geophysics  Directorate,  from  1.2  to  30  km, 
Holloman  Air  Force  Base,  NM,  7  December  1990  (Probe  Set  2) . 


ALT 

<Km)  1 5- 


7igur«  €.  Log  (Raw  c/}  vs.  Altitude  Profile  (Night)/ 
Collected  by  the  Geophysics  Directorate,  from  1.2  to  30  km, 
Holloman  Air  Force  Base,  NM,  8  December  1990  (Probe  Set  1} . 


24 


rlgur*  7. Log (Raw C/j vs.  Altitude  Profile  (Night) , 
Collected  by  the  Geophysics  Directorate,  from  1.2  to  30  kn, 
Holloman  Air  Force  Base,  NM,  8  December  1990  (Probe  Set  2) . 


b.  7  D9CBmb9r  1990 


On  7  December  1990,  a  day  flight  with  both  parallel  and 
perpendicular  probe  configurations,  the  two  different  sets  of  probes  behaved 
differently.  The  set  of  parallel  probes  exhibited  a  trend  common  to  the  other  day 
flights.  At  low  altitudes,  up  to  approximately  8-12  km,  some  of  the  data  points 
closely  approached  the  baseline  temperature-pressure  cun/e.  As  the  probe 
continued  to  ascend,  the  log  (Raw  C„^  values  Increased  and  moved  away  from 
the  system  noise  level  curve  by  an  order  of  magnitude  (Figure  4).  The 
perpendicular  set  of  probes  exhibited  enhanced  differences.  These  log  (Raw  C/) 
values  consistently  were  an  order  of  magnitude  higher  than  the  noise  level 
throughout  the  entire  flight  (Figure  5).  This  indicates  that  there  was  a  constant 
source  of  error,  or  contamination,  throughout  the  entire  flight. 

Rgures  8  thru  1 1  show  the  x-y  correlation  between  the  two  values 
(CT1  and  CT2)  measured  at  each  altitude.  These  figures  show  the  effect  of  noise 
or  some  constant  signal  with  a  value  of  approximately  0.006  K  that  affects  the 
data  in  channel  2.  Examination  of  10  km  increments  of  the  atmosphere  reveals 
that  this  signal  first  became  apparent  between  10-20  km  (Figure  10)  and  became 
very  pronounced  at  20-30  km  (Rgure  11).  The  dense  spike  in  the  data  at 
CT1  *0.006  K  Illustrates  this  effect. 
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Figura  8.  CTl  vs.  CT2  (Day),  correlation  of  RMS  Temperature 
Fluctuations  Between  Two  Pairs  of  Probes 
(Parallel/Perpendicular  Configuration)  Collected  Between  1.2 
and  30  km,  Holloman  Air  Force  Base,  NN,  7  December  1990. 
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Vigvatm  9.  CTl  vs.  CT2  (Day),  Correlation  of  RMS  Temperature 
Fluctuations  Between  Two  Pairs  of  Probes 
(Parallel/Perpendicular  Configuration)  Collected  Between  1.2 
and  10  km,  Holloman  Air  Force  Base,  NM,  7  December  1990. 
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rigur*  lO.  CTl  vs.  CT2  (Day),  correlation  of  RMS  Temperature 
Fluctuations  Between  Two  Pairs  of  Probes 
(Parallel/Perpendicular  Configuration)  Collected  Between  10 
and  20  km,  Holloman  Air  Force  Base,  NM,  7  December  1990. 
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Figiurs  11.  CTl  vs.  CT2  (Day) ,  Correla'tion  of  RMS  Temperature 
Fluctuations  Between  Two  Pairs  of  Probes 
(Parallel/Perpendicular  Configuration)  Collected  Between  20 
and  30  km,  Holloman  Air  Force  Base,  NM,  7  December  1990. 


The  correlation  coefficient  for  the  entire  flight  was  50%,  the  lowest 
of  the  five  successful  flights.  Table  II  shows  that  r  was  high,  96%,  in  the  first  5 
km  interval  and  was  below  50%  for  the  rest  of  the  flight. 


Table  11.  CT1  vs.  CT2  Correlation  for  7  December  1990  Daytime  Launch  with 
Parallel/Perpendicular  Probe  Configuration  at  Holloman  Air  Force  Base,  NM. 


ALTITUDE  (km) 

CORRELATION,  r 

Data  points 

1.2-5 

.96 

155 

5-10 

.42 

220 

10-15 

.32 

252 

15-20 

.27 

174 

20-25 

.34 

240 

25-30 

.40 

176 

TOTAL 

.50 

1217 

c.  8  December  1990 

The  plots  of  log  (Raw  C„*)  values  for  the  8  December  1990  flight 
followed  the  characteristic  nighttime  behavior.  Some  of  the  data  consistently 
approached  the  noise  level  up  to  approjdmately  15-20  km,  at  which  point  they 
showed  a  slight  shift  to  the  higher  values.  Figures  6  and  7  show  these  results. 


31 


Figures  12  thru  15  plot  the  CT1  vs.  CT2  correlation  vs.  altitude. 
These  figures  show  the  effect  of  a  small  but  relatively  constant  offset.  The  CT2 
vedues  are  consistently  greater  than  the  CT1  values.  This  could  not  occur  without 
a  signal,  instrument,  or  csdibration  error  since  the  probes  passed  through 
essentially  the  same  air  mass  up  to  30  km  and  took  893  independent 
measurements.  Tsdsle  III  shows  that  the  correlation  coefficient  was  high  in  ail 
intervals  with  the  exception  of  the  5-10  km  region  which  had  an  insufficient 
number  of  data  points  above  the  noise. 


TABLE  III.  CT1  vs.  CT2  Correlation  for  Geophysics  Directorate  8  December  1990 
Night  Launch  with  Dual-Parallel  Probe  Configuration  at  Holloman  Air  Force  Base, 
NM. 


ALTITUDE  (km) 

CORREUTION,  r 

Data  Points 

1.2-5 

.98 

115 

5-10 

.58 

29  * 

10-15 

.90 

161  1 

15-20 

.92 

240  1 

20-25 

.86 

199  1 

25-30 

.86 

148 

TOTAL 

.91 

893 

*  statistically  insufficient 
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Rgur*  12.  CT1  vs.  CT2  (Night),  Correlation  of  RMS  Temperatijre  Ructuations 
Between  Two  Pairs  of  Probes  (Dual-Parallel  Configuration)  Collected  Between 
1.2  and  30  km.  Holloman  Air  Force  Base,  NM,  8  December  1990. 
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Rgur»  13.  CT1  vs.  CT2  (Night),  Correlation  of  RMS  Temperature  Ructuations 
Between  Two  Pairs  of  Probes  (Dual-Parallel  Configuration)  Coiiected  Between 
1.2  and  10  km,  Holloman  Air  Force  B€»e.  NM,  8  December  1990. 
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Rgitr*  14.  CT1  vs.  CT2  (Night),  Correlation  of  RMS  Temperature  Ructualions 
Between  Two  Pairs  of  Probes  (Dual-Paraiiel  Configuration)  Collected  Between 
10  and  20  km,  Holloman  Air  Force  Base,  NM,  8  December  1990. 
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Between  Two  Pairs  of  Probes  (Dual-Parallel  Configuration}  Collected  Between 


20  and  30  km,  Holloman  Air  Force  Base,  NM,  8  December  1^. 


36 


The  plots  of  the  4,  5,  and  6  December  flights  followed  the  typical 
daytime  trend  and  are  in  Appendix  A.  The  data  collection  during  the  descent 
phase  of  the  experiment  was  not  successful  due  to  difficulties  in  using  the  heavier 
than  standard  package. 

4.  Analysis 

The  high  degree  of  correlation  of  ail  flights  in  the  1-5  km  region  is  due 
in  part  to  the  turbulent  boundary  layer  which  produces  high  in  the  daytime. 
Following  Weitekamp,  solar  induced  temperature  differences  due  to  insolation  are 
smallest  at  ground  level  and  are  the  largest  at  the  tropopause  when  the  probe 
is  at  a  temperature  of  approximately  200  K.  During  the  day,  solar  energy  raises 
the  probe  surface  temperature  above  the  ambient  air  temperature.  This 
difference  increases  as  the  probe’s  temperature  declines.  This  increase  mimics 
turbulent  fluctuations  when  it  varies  rapidly  with  time.  At  night,  however,  the 
surface  temperature  drops  below  that  of  the  air  as  it  cools  by  radiation.  The 
decrease  in  correlation  in  the  5-10  km  region  appears  to  be  the  result  of  the 
decline  in  the  boundary  layer  contribution  coupled  with  the  increase  in  the  solar 
induced  probe-atmosphere  temperature  differences. 

The  only  flight  with  the  parallel/perpendicular  configuration  of  the  two 
sets  of  probes,  launched  on  7  December,  had  by  far  the  lowest  correlation  (total 
average  r  of  50  %).  This  indicates  that  the  probe  orientation  with  respect  to  the 
sun  influenced  the  differences  in  the  measured  temperatures  [Ref.  3].  The 
motion  of  the  ascending  sonde  package  continuously  changed  the  orientation. 
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and  hence  the  projected  surface  area,  of  the  probes  respect  to  the  incident 
solar  radiation  resulting  In  uneven  thermodynamic  heating  and  cooling  of  the 
probes  (Ref.  3;p.  33].  The  results  of  the  7  December  1990  flight  further  support 
his  findings. 

The  most  significant  results  occurred  in  the  region  of  most  pressing 
Interest,  namely  10  to  20  km.  The  nighttime  data  correlation  was  approximately 
91  %,  markedly  higher  than  any  of  the  daytime  data  (Rgures  10  and  14).  This 
indicated  that  solar  contamination  of  the  daytime  data  is  a  real  effect. 

The  4  December  1990  flight  had  very  low  correlation  in  the  15  km  to 
25  km  region;  much  lower  than  even  the  7  December  flight  (see  Table  IV).  This 
could  have  been  due  to  instrument  error,  calibration,  or  perhaps  the  probes 
actually  measured  an  increase  in  turbulence.  Some  sort  of  instrument 
malfunction  was  most  likely  the  case.  The  fact  that  the  correlation  of  this  flight 
jumped  to  77%  from  25  km  to  30  km,  compared  to  16%  in  the  20  to  25  km 
altitude  region  indicates  a  systematic  error.  This  marked  increase  was  the  most 
dramatic  data  discontinuity  between  any  altitude  regions  in  any  of  the  five  flights 
in  Diurnal  Variation  Experiment  2  and  should  be  considered  skeptically  prior  to 
making  any  inferences  from  it. 

Rgure  16  is  a  plot  of  the  results  of  the  correlation  analysis  of  each 
flight.  T^le  IV  summarizes  the  data  for  the  4,  5,  and  6  December  flights. 
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Rgura  16.  Correlation  of  CT1  and  CT2  vs.  Altitude  for  ail  Diurnal  Variation 
Experiment  Rights  Conducted  by  the  Geophysics  Directorate  at  Holloman,  Air 
Force  Base,  NM,  from  4-8  December  1990. 
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TABLE  IV.  CT1  and  CT2  Correlation  vs.  Altitude  for  Probe  Pairs  for  Remaining 
Rights  in  Diurnal  Variation  Experiment  2. 
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C.  OVE  SUMMARY 


The  high  correlation  obtained  for  the  temperature  differences  (CT1  and  CT2) 
for  the  sets  of  probes  in  the  night-time  flight  on  8  December  1990  and  the  low 
con'elation  obtained  during  day-time  flights  particularly  near  the  tropopause 
support  the  hypothesis  that  the  apparent  order  of  magnitude  increase  in  optical 
turbulence  is  an  instrument  effect  attributable  to  solar  heating  of  the  probes,  lliis 
is  consistent  with  previous  work  that  shows  that  the  solar  induced,  probe-air 
temperature  differences  are  largest  at  low  temperatures  [Ref.  3:p.  19].  The  fact 
that  the  only  flight  that  had  the  parallei/perpendicular  configuration  of  the  two 
sets  of  probes,  the  7  December  launch,  had  the  largest  overall  variations 
(smallest  correlation  coeffident,  r),  indicsttes  that  the  probe  orientation  with 
respect  to  the  sun  contributed  to  the  differences  in  the  measured  temperatures. 

The  motion  of  the  package  (rotation  or  pendulum-like  osdllations),  could 
produce  a  low  frequency  signal.  The  parallei/perpendicular  configuration  would 
correlate  more  poorly  than  the  double  parallel  configuration  because  the  probes 
would  move  through  the  air  masses  with  different  orientations  and  be  exposed 
to  different  conditions.  Differential  solar  heating  of  the  probes  from  probe  vwre 
oscillation,  possibly  exacerbated  by  payload  motion  (pendulum  and  rotation) 
could  produce  a  higher  frequency  signal  which  could  mask  the  data  [Ref.  3:pp. 
16,46].  The  flexible  probe  wires,  which  oscillate  from  the  cross  flow  of  air,  would 
continuously  present  differing  cross-sectional  areas  to  the  sun  and  induce  a 
signal.  If  probe  vibrations  were  the  cause  of  these  differing  results,  then  both  the 
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parallel/perpendicular  configuration  and  the  dual  parallel  configuration  would 
have  exhibited  uneven  heating  behavior.  The  daytime  data  currently  available  are 
only  reliable,  meaning  consistent,  with  nighttime  data  up  to  approximately  8-10 
km,  which  is  an  cdtitude  that  is  clearly  beneath  the  tropopause.  Daytime 
microthermal  probe  data  beyond  this  are  highly  questionable.  Although  the 
number  of  flights  in  this  series  of  launches  was  limited,  these  results  are 
consistent  with  those  of  many  previous  day  launches  and  support  the  hypothesis 
that  solar  heating  of  the  probes  contaminates  the  daytime  data.  Although  the 
exact  mechanism  of  this  contamination  is  not  readily  apparent,  the  evidence 
confirming  its  existence  is  convincing. 

if  balloon  wake  (thermal  shedding)  was  responsible  for  the  diurnal  variation 
in  the  C/  measurements,  the  correlation  of  the  nighttime  flights  should  be  poorer 
than  the  day.  The  balloon,  which  is  cooler  than  its  surroundings  at  night,  creates 
a  cooled  wake  beneath  it  as  it  rises.  The  wake  effect  should  be  worse  at  night 
since  the  temperature  difference  between  the  balloon  and  the  surrounding  air  is 
larger  than  in  the  day.  Weitekamp  showed  that  this  effect  was  small  at  flie 
distances  that  the  thermosonde  packages  were  suspended  beneath  the  balloon 
[Ref.  3:p.  25]. 

D.  SONOE  STABiUZATiON  FOR  SUBSEQUENT  EXPERIMENTS 

To  eliminate  solar  effects  apparent  in  the  daytime  leujnches  due  to  rotation 
and  oscillation  of  the  sonde,  it  is  necessary  to  stabilize  the  sonde  ^  much  as 
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possible.  The  most  reliable  and  consistent  data  will  occur  when  the  probes 
maintain  a  consistent  alignment  with  respect  to  the  sun.  Possible  solutions  to 
this  may  include  the  use  of  a  gyro  mechanism  to  stabilize  the  platform,  vertical 
baffles  or  fin  stabilizers,  solar  cells  to  determine  relative  solar  angle,  and 
weights/anchors  to  hold  the  apparatus  steady.  To  reduce  the  effects  of  probe 
vibration,  the  tungsten  wire  probes  must  be  attached  straight  across  the  probe 
supports  rather  than  with  a  flexible  configuration.  This  will  involve  a  delicate 
compromise  to  avoid  thermal  stress  and  strain  on  the  wire  and  its  contacts  and 
yet  suppress  vibration. 
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iV.  PROBE  RESPONSE 


A.  PROBE  RESPONSE  TIME  AS  A  FUNCTION  OF  DECREASING 
TEMPERATURE  AND  AIR  DENSITY 

1.  Th»nnocoupl«  MMSur«m«nt  Limitations  Dus  to  Probs  Rssponsa 
Tima 

A  second  systematic  issue  that  Impacts  the  performance  of  the  thermal 
probe  turbulence  measurements  involves  estimating  the  probe  response  time  vs. 
altitude.  When  a  thermocouple  measures  the  temperature  of  its  atmospheric 
surroundings,  the  value  that  it  provides  is  its  own  temperature,  not  that  of  the 
surrounding  air.  This  occurs  because  of  the  finite  heal  transfer  mechanisms 
occurring  at  the  surface-air  interface.  The  themnaJ  mass  of  the  probe  and  the 
convective  air-surface  interface  limit  the  rate  of  heat  transfer  producing  a  delayed 
probe  response.  This  limits  the  high  frequency  components  of  the  signal  used 
to  compute  structure  functions.  [Ref.  15;p.  15]  This  becomes  increasingly 
important  for  structure  function  measurements  with  small  separation  distances. 
A  10  cm  separation  requires  a  probe  response  time  of  1-2  ms.  This  problem 
increases  at  higher  altitudes.  The  low  density  environment  in  the  1 0  to  20  km  (1 0 
mb)  region  offers  less  effective  heat  transport  (primarily  forced  convection)  by 
which  the  probes  adjust  to  the  temperature  of  their  surroundings,  consequently 
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measurements  underestimate  the  corresponding  structure  functions  as  altitudes 
increase. 

2.  Fluid  Characteristics  and  Heat  Transfer 
a.  Forced  Convection 

Forced  convection  is  the  dominant  (albeit,  still  a  weak)  heat 
transfer  mechanism  in  regions  where  the  mean  free  path  of  the  surrounding 
molecules  is  comparable  to  the  dimensions  of  the  heat  transfer  surface  (the 
cylindrical  wire  probe).  The  physical  mechanism  of  forced  convection  which 
removes  heat  (energy)  from  the  probes  is  the  collisional  momentum  transfer  to 
surrounding  air  molecules  when  these  molecules  impact  the  probe  [Ref.  16:p. 
6311.  mean  free  path  of  atmospheric  molecules  increases  from  0.06  fjm  at 
sea  level  to  approximately  4.4  //m  at  30  km  [Ref.  17:p.  14-13].  This  is 
compar^le  to  the  dimensions  of  a  12.4  /vm  diameter  cylindrical  wire  probe. 
These  are  the  dimensions  of  the  ANSI  Type  E  chromel-constantan  thermocouples 
that  the  NFS  group  used  in  the  summer  and  fall  1992  experiments:  the 
thermocouples  were  manufactured  by  Omega  Engineering,  Inc.  [Ref.  18:p.  A-33]. 
Above  30  km,  as  the  moleculau*  mean  free  path  length  continues  to  increase, 
forced  convection  is  no  longer  the  dominant  heat  transfer  mechanism  for  probes 
of  this  size. 
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b.  WscosAy 

The  thermodynamic  transport  property  that  characterizes  the  fluid- 
mechanical  behavior  of  the  medium  is  the  viscosity,  //,  aJso  known  as  the 
dynamic  viscosity.  It  exists  as  a  point  function  in  a  flowing  fluid,  the  air,  and 
follows  all  of  the  laws  and  state  relations  of  ordinary  equilibrium  thermodynamics 
[Ref.  7:p.  20].  Viscosity  is  a  measure  of  the  diffusion  of  momentum  parallel  to  the 
flow  velocity  and  transverse  to  the  gradient  of  the  flow  velocity  [Ref.  16:p.  219] 
and  has  the  dimensions  of  (kglm-s).  It  relates  the  local  stresses  in  a  moving 
fluid  to  the  rate  of  strain  in  the  fluid  element  by 


du 

dy' 


(13) 


where  r  is  the  shear  stress  and  du/dy  is  the  local  fluid  velocity  gradient  [Ref.  7;p. 
26], 

Prandtl  showed  that  low  viscosity  fluid  flows  such  as  airflow  are 
divided  into  a  thin  viscous  boundary  layer  near  a  solid  surface  which  borders  a 
nearly  inviscid  outer  layer  which  is  still  governed  by  the  Euler  and  Bernoulli 
Equations  [Ref.  19:p.  46]  Under  shearing,  a  fluid  begins  to  move  at  a  strain  rate 
inversely  proportional  to  the  coefficient  of  viscosity.  Viscosity  depends  primarily 
on  temperature  at  low  altitudes  by 


P.2’2/3 
T*S  ' 


(14) 
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where  is  a  constant  equsy  to  1.458  x  lO"*  kgs’m  ’IC’'®  and  S,  Sutherland’s 
constant,  is  equal  to  110.4  K  [Ref.  20:p.  14-5]. 

c.  KiMmttIc  VI»co9tty 

The  kinematic  viscosity,  v,  is  the  ratio  of  the  dynamic  viscosity  of 
a  gas  to  its  density  by 


(15) 


where  p  is  the  atmospheric  density  in  (kg/m*)  [Ref.  20;p.  14-6].  The  dimensions 
of  V  are  (m*/s).  Figure  17  shows  how  kinematic  viscosity  increases  with  altitude. 


KINEMATIC  VISCOSITY  VS.  ALT 


Rgure  17.  Kinematic  Viscosity  as  a  Function  of  Altitude  [Ref.  20:p.  14-6}. 


ef.  No^llp  Condition 

Vy^en  a  fluid  flow  encounters  a  solid  surfsu^e,  molecular 
interactions  cause  the  fluid  and  the  surface  to  seek  energy  md  momentum 
equilibrium  with  each  other;  all  liquids  are  essentially  in  equilibrium  vsnth  the 
surface  with  which  they  are  in  contact;  fluids  in  contact  with  a  solid  take  on  the 
velocity  and  temperature  of  the  surface.  These  are  the  no-siip  and  no¬ 
temperature-jump  conditions  which  are  only  valid  when  the  molecular  mean  free 
path  is  less  than  the  surface  dimensions  of  the  solid.  [Ref.  7;p.  35] 

e.  Reynold*  Number,  Re 

The  parameter  that  describes  the  viscous  behavior  of  a  Newtonian 
fluid  such  as  air  is  the  dimensionless  Reynolds  number  Re 

(16) 

|A  V 

where  p  is  the  atmospheric  density,  p  is  the  dynamic  viscosity,  V  is  the 
characteristic  velocity,  and  L  is  the  characteristic  length  scale  of  the  flow.  The 
Reynolds  number  is  a  ratio  of  the  dynamic  forces  of  mass  flow  (inertia)  to  the 
viscous  forces  within  the  atmosphere.  For  a  12.4  pm  chromel-constantan  probe, 
the  value  of  the  Reynolds  number  varies  from  3.75  at  sea  level  to  0.07  at  30  km 
for  vwndspeeds  of  4.4  m/s.  Figure  18  is  a  plot  of  the  Reynolds  number  for  the 
12.4  pm  probe.  The  plot  includes  Reynolds  numbers  computed  using  both  the 
maximum  and  minimum  probable  diameters  of  a  12.4  pm  wire  in  accordance  with 
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the  manufacturer’s  manufacturing  tolerances  (a  12.4  //m  AWG  No.  56 
thermocouple  wire  can  vary  from  11.4  to  13.4  //m  [Ref.  21]).  These  values 
represent  two  regimes  of  the  Reynolds  number;  0<Re<1.0  describes  a  region 
of  highly  viscous  laminar  creeping  motion,  and  the  range  1 .0<Re<100  describes 
a  region  of  laminar  motion  with  strong  Reynolds  number  dependence  [Ref.  7:  p. 
309],  Table  V  summarizes  Reynolds  numbers  for  the  probe  as  a  function  of 
altitude  where  the  probe  velocity  is  4.4  m/s. 


TABLE  V.  Reynolds  Number  for  NFS  12.4  +/-  1,0  //m  Chromel-Constantan 
Probes  from  0  to  30  km  at  v=4.4  m/s. 


ALT  (km) 

1 1 .4  //m  Probe 

12.4  //m  Probe 

13,4  //m  Probe  | 

0 

3.44 

3.75 

4.05  1 

10 

1.43 

1.55 

1.68  1 

20* 

0.32 

0.34 

0.37  1 

30 

0.06 

0.07 

0.07  1 

*  above  19  km,  where  Re  <0.4,  the  equation  for  Re  begins  to  lose  its 
accuracy;  however,  empirical  fits  to  the  data  by  Holman  [Ref.  16:p.  303] 
have  shown  that  Re  values  follow  approximately  the  same  trend  in  this 
regime. 
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REYNOLDS  NUMBER  (Re) 


REYNOLDS  NUMBER  VS.  ALTITUDE 

12.4  Micron  Probe  {V=4.4  m/s) 


Rgur«  18:  Reynolds  Number  vs.  Altitude  for  NFS  12.4  +/-  1.0  //m  Chromel- 
Constantan  Thermocouple  Probe. 
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f.  Nuss9lt  Number 


From  Kothandaraman  and  Subramanyan  the  Nusselt  number  for 
a  cylinder  in  a  cross-flow  is  [Ref.  22:p.  91} 


where  C  and  n  are  empirically  developed  constants  based  upon  the  local 
Reynolds  number  where 


(18) 


Table  VI  contains  the  values  of  the  constants  C  and  n.  Rgure  19  is  a  plot  of  Nu 
vs.  altitude  for  the  NFS  probe  based  on  the  Reynolds  numbers  shown  in  Table 
V.  Tedsle  VII  summarizes  the  Nusselt  number  values. 
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TABLE  Vi.  Constants  for  Determination  of  Nusselt 
Number  for  Eq.  (16). 


Re 

C 

n 

<  0.4 

0.891* 

0.330* 

0.4  -4.0 

0.891 

0.330 

>  4.0 

0.821 

0.385 

*  empirically  fit 


TABLE  Vli.  Nusselt  Number  vs.  Altitude  for  NPS  12.4  +/- 1.0  //m  Chromel- 
Constantan  Thermocouple  Probes. 


ALT  (km) 

1 1 .4  ^m  Probe 

12.4  nm  Probe 

13.4  //m  1 

Probe  1 

0 

1.34 

1.38 

1.41 

10 

1.00 

1.06 

20 

0.61 

0.63 

0.64 

30 

0.35 

0.37 

0.38 
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NUSSELT  NUMBER 


NUSSELT  NUMBER  VS.  ALTITUDE 

12.4  Micron  Probe  (V=4.4  m/s) 


1.6 
1.4 
1.2 
1 

0.8 
0.6 
0.4 
0.2 
0 

0  5  10  15  20  25  30 

ALTITUDE  (km) 


Rgura  19.  Nusselt  Number  vs.  Altitude  for  NFS  12.4  +/-  1.0  Chromel- 
Constantan  Thenmocouple  Probe. 
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g.  ffat«  of  ConvctIVB  Hoat  Flow 


Newton's  Law  of  Cooling  governs  the  rate  of  convective  heat  flow 
between  a  solid  and  its  environment  by  [Ref.  19:p.  89] 

q^M^T,  (19) 


where 

q  -  rate  of  heat  transfer 
h  =  convective  heat  transfer  coefficient 
A  =  solid  surface  area 

AT  =  temperature  difference  between  the  solid  and  its  environment. 

h.  ConvacUva  Haat  Tranafar  Coafficlant, 

The  convective  heat  transfer  coefficient,  h^,  is  a  parameter  that 
describes  the  energy  transport  (both  thermodynamic  and  mass)  between  the 
probe  and  the  surrounding  atmosphere,  it  is  a  fluid  property  that  conveniently 
groups  terms  that  are  a  complex  function  of  fluid  flow  and  geometric  conditions 
of  flie  probe  and  the  environment  [Ref.  19:p.  44].  It  is  generally  expressed  as  an 
empirical  form  of  the  Nusselt  number,  Nu,  by 


where  d  is  the  characteristic  length  (wire  diameter)  and  k  is  the  thermal 
conductivity  of  the  surrounding  air  in  (kcal/m-hr  K)  [Ref.  23:p.  6-3].  The 
convective  heat  transfer  coefficient,  h^,  has  dimensions  of  W/m*K, 

Previous  work  by  Brown  [Ref  12:p.  19]  and  Weitekamp  [Ref.  3:p. 
38]  has  shown  that  there  is  approximately  a  factor  of  four  decrease  in  the 
convective  heat  transfer  coefficient,  h^,  from  the  ground  to  an  altitude  of  30  km. 
The  value  of  h^  for  the  NPS  12.4  fjm  chromel-constantan  probe  exhibited  similar 
behavior;  the  actual  values  of  h^  were  different  due  to  the  different  sizes  and 
composition  of  the  respective  probes,  but  the  trend  was  the  same  as  shown  in 
Figure  20  for  the  NPS  probes  ascending  at  4.4  m/s.  Table  VIII  summsuizes  these 
data 


TABLE  VIII.  Convective  Heat  Transfer  Coefficient,  h^,  vs.  Altitude  for  NPS  12.4 
+/- 1.0 //m  Chromel-Constantan  Thermocouple  Probes. 


ALT  (km) 

h^  for  1 1 .4  //m 
Probe  (W/m®K) 

h,.  for  12.4//m 
Probe  (W/m*K) 

he  for  13.4  //m 

Probe  (W/m*K) 

0 

2973 

2808 

2665 

10 

1760 

1662 

1577 

20 

1040 

982 

932 

30 

637 

602 

571 
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4000 


HEAT  TRANSFER  COEFFICIENT,  h 

12,4  Micron  Probe  (V=4.4  m/s) 


3500 


3000 


Flgur«  20.  Atmospheric  Convective  Heat  Transfer  Coefficient,  h^,  vs.  Altitude  of 
NPS  12.4  +/-  1.0  /mi  Chromel-Constantan  Thermocouple  Probe  with  VMnd 
Velocity  of  4.4  m/s. 
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B.  TIME  CONSTANT,  r 


1.  Theoretical  Probe  Response 

The  time  for  the  thermocouple  to  complete  63.2  %  of  its  response  to 
a  step  change  in  the  temperature  of  the  air  around  it  is  the  characteristic  time,  r 
[Ref.  24:  p.  566].  The  characteristic  time  is 

(21) 


where  p  is  the  thermocouple  density  In  (gm/cm®),  c  is  the  thermocouple  specific 
heat  capacity  in  (J/g*K),  d  is  the  diameter  of  the  thermocouple  in  (cm),  and 
is  the  convective  heat  transfer  coefficient  in  {J*s)/(cm**i<)  Hence,  the 
characteristic  time  is  the  ratio  of  the  probe  thermal  storage  capacity  to  the  heat 
input  per  degree  of  temperature  difference.  (Ref.  24:p.  555] 

Eq.  (21)  shows  that  the  response  time  of  a  particular  probe  varies  only 
with  the  convective  heat  transfer  coefficient,  h^;  the  other  quantities  depend  on 
the  structure  of  the  thermocouple  and  do  not  change  appreciably  under 
experimental  conditions.  Several  factors  affect  the  accuracy  of  this  equation. 
Manufacturing  tolerances  in  wire  diameter  coupled  with  the  presence  of  a  welded 
junction  bead  that  connects  the  thermocouple  wires  near  the  center  (Figure  21) 
contribute  to  an  imprecise  prediction  of  probe  response. 
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IDEAL  WELDED 
JUNCTION  BEAD 


D  =  3xd 


Rgura  21.  Thermocouple  Junction  Bead. 
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The  geometry  and  mass  of  this  welded  junction  bead  give  it  different 
heat  transfer  characteristics  than  those  of  the  cylindrical  wire,  and  therefore,  affect 
the  true  value  of  t.  Moffat  accounted  for  this  difference  by  the  expression 


T  (Df^ 


(22) 


where  r  is  the  true  characteristic  time,  Tq  is  the  characteristic  time  when  the 
diameter  of  the  junction  bead  and  the  diameter  of  the  wire  are  the  same  (no  weld 
effect),  0  is  the  diameter  of  the  bead  (assumed  to  be  spherical),  and  d  is  the 
diameter  of  the  probe  wire.  Moffat  also  emphasized  that  to  ensure  first  order 
probe  response,  the  diameter  of  the  weld  should  not  exceed  the  wire  diameter 
by  more  than  10%.  [Ref.  24:p.  566] 

2.  Thermocouple  Characteristics 

(1)  Thermocouple  Type 

The  NFS  differential  thermosonde  used  thermocouple  rather 
than  resistance  probes  to  simplify  the  electronic  circuit.  The  smallest 
commercially  available  thermocouple  was  12,4  ^m  in  diameter,  which  was  four 
times  larger  than  those  used  by  the  Geophysics  Directorate.  The  increased 
thermal  mass  reduced  the  high  frequency  response  which  had  increasing 
importance  in  structure  function  measurements  for  small  separations  ( ~  10  cm). 
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(2)  Material  Properties 


The  NPS  Group  used  ar>  ANSJ  Type  E  cylindrical  wire 
thermocouple  to  make  temperature  measurements  for  all  launches  from  July  to 
November  1992.  According  to  the  specifications  of  the  manufacturer,  Omega 
Engineering,  Inc.,  this  chromel-constantan  thermocouple  was  12.4  pm  '.n 
diameter.  Its  composition  was  75%  (Ni-10Cr)  and  25%  (Cu-43Ni)  which  gave  it 
a  specific  heat  capacity,  c,  of  .4367  (J/g  •  K).  Appendix  C  has  more  details  on  the 
thermocouple  composition  and  characteristics.  The  NPS  Group  chose  this 
thermocouple  since  it  was  the  smallest  available  and  because  of  its  large 
thermoelectric  response. 

(3)  Thermocouple  Geometry 

To  ensure  that  all  consteaits  affecting  the  time  constant  were 
valid,  we  measured  the  diameter  of  the  chromel-constantan  thermocouple  wire 
with  a  Cambridge  200,  tungsten  filament  scanning  electron  microscope  (SEM). 
The  measured  diameter  of  the  cylindrical  segment  of  the  wire  was  13.7  prr\  with 
an  unspecified  tolerance.  This  value  was  10%  higher  than  the  manufacturer's 
specifications  for  the  AWG  no.  56  thermocouple  wire.  Since  the  SEM  had  not 
been  calibrated  recently  and  the  operator’s  manual  contained  no  specific 
measurement  tolerances,  we  used  the  measurement  as  a  verification  of  the 
general  dimensions  of  the  wire  and  to  estimate  the  bead  size  [Ref.  25]. 

The  V-shap3d  wire  had  a  welded  bead  near  the  v,  which 
according  to  the  manufacturer,  should  have  been  spherically  shaped  and 
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approximately  3  times  the  wire  diameter  [Ref.  18:p.  A-31}.  The  SEM  measured 
the  dimensions  of  the  roughly  cylindrical  bead  to  be  approximately  15.5  //m  x 
15.5  tjn\  X  52.3  //m  resulting  in  a  volume  of  approximately  13, OCX)  /ym^.  Since  the 
SEM  only  measured  in  2  dimensions,  it  did  not  provide  a  value  for  the  bead's 
height.  Photograph's  of  the  bead  showed  that  the  center  segment  of  the  bead 
appeared  to  be  cylindrical;  therefore,  we  assigned  the  bead  a  thickness  of 
15.5/;m,  the  same  value  as  the  width.  We  used  a  spherical  approximation 
because,  ideally,  this  would  be  the  shape  of  the  weld  bead  according  to 
manufacturer's  specifications  and  as  found  in  the  examination  of  numerous  other 
thermocouple  wires.  By  approximating  the  shape  of  the  bead  as  spherical  and 
using  the  standard  equation  for  the  volume  of  a  sphere, 

(23) 

6 

where  D  is  the  diameter  of  the  sphere,  the  effective  diameter  of  the  bead  became 
29  fjTT).  This  diameter  was  -2.3  times  the  wire  diameter,  in  reasonable 
agreement  with  the  3:1  ratio  given  by  the  manufacturer.  We  calculated  an 
adjusted  value  for  the  bead  diameter  due  to  the  discrepancy  between  the 
manufacturer’s  specifications  and  the  measured  SEM  vaUues  for  the  wire 
diameter.  Since  the  SEM  value  was  10%  higher  than  the  manufacturer’s  value 
for  diameter,  its  values  for  the  bead  dimensions  should  have  been  10%  higher 
as  well.  Factoring  this  adjustment  into  Eq.  (23),  yielded  an  adjusted  bead 
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diameter  of  26  //m.  Using  either  the  SEM  or  manufacturer’s  values,  the  ratio  of 
the  bead  diameter  to  the  wire  diameter  dearly  exceeded  the  1 0%  allowance  for 
first-order  effects,  therefore,  higher-order  effects  would  affect  the  probe  response. 

3.  Theoretical  Calculation  of  r 

Figure  22  shows  a  plot  of  r  vs.  altitude  for  the  NPS  12.4  ^m  diameter 
wire  thermocouple,  uncorrected  for  a  junction  bead.  The  uncorrected  r  for  the 
12.4  ywm  probe  had  a  value  of  4.2  ms  at  0  km  and  19.7  ms  at  30  km.  Figure  23 
shows  a  plot  of  T  vs.  altitude  for  the  NPS  12.4  //m  diameter  wire  thermocouple, 
corrected  for  a  26  //m  junction  bead.  The  time  constant  accounting  for  the  bead 
was  5.8  ms  at  0  km  and  27.0  ms  at  30  km.  The  effect  of  correcting  for  the  bead 
yielded  response  times  that  were  ~30%  slower  than  for  the  calculations  with  no 
bead  correction.  Table  iX  summarizes  these  data. 


TABLE  IX.  Calculated  Thermocouple  Probe  Response  vs.  Altitude,  for  a  12.4  +/- 
1.0  //m  Chromel-Constantan  Wire,  with  and  without  Correction  for  a  26  //m 
Junction  Bead. 


d  ijjm) 

11.4 

11.4 

12.4 

12.4 

13.4 

13.4  1 

ALT 

(km) 

no 

bead 

26  ^m 
bead 

no  bead 

26  //m 
bead 

no  bead 

1 

26  ^m 
bead 

0 

EBH 

o 

in 

4.2 

5.6 

4.8 

6.2 

10 

6.2 

8.4 

7.1 

9.4 

8.1 

10.4 

20 

10.5 

14.3 

12.1 

15.9 

13.8 

17.6  1 

30 

17.1 

23.3 

19.7 

26.0 

22.5 

28.8  1 

62 


TIME  CONSTANT  (s) 


Hgure  22.  Calculated  Time  Constant  for  NPS  12.4  +/- 1.0  fjm  Type  E  Chromel> 
Constantan  Thermocouple  Probe  with  no  Correction  for  a  Junction  Bead. 
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Constantan  Thermocouple  Probe  with  a  Correction  for  a  26  //m  Junction  Bead. 
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C.  EXPERIMENTAL  VERIFICATION  OF  PROBE  RESPONSE  TIME 


1.  Empirical  Formula  For  Probe  Thermal  Response 

When  used  to  measure  local  gas  temperature,  a  bare-wire 
thermocouple  displays  physical  and  thermodynamic  behavior  similar  to  that  of 
a  cold-wire  anemometer.  A  cold-wire  anemometer  measures  rapidly  fluctuating 
flows  such  as  the  turbulent  boundary  layer  surrounding  an  object  (the  probe). 
It  consists  of  a  fine  wire  mounted  between  two  small  probes.  An  empirical 
equation  that  relates  the  wire  response  time  with  the  windspeed  is 

(24) 

X 

where  a  and  b  are  constants,  a  being  the  y-intercept  of  the  plotted  data  and  b 
being  the  slope  of  this  data,  plotted  according  to  linear  regression  where  r  is  the 
response,  and  v  is  tfie  wind  velocity  (Ref.  26:p.  13]. 

2.  Wind  Tunnel  Test 

We  constructed  a  wind  tunnel  and  subjected  a  12.4  //m  probe  (not  the 
same  probe  that  we  had  SEM  diameter  measurements  on)  to  cross-winds  varying 
in  speed  from  0  to  4.4  meters-per-second.  The  maximum  vwndspeed  was  4.4 
m/s  because  of  the  operating  limit  of  the  24V  DC  fan  that  was  the  wind  source. 
The  4.4  m/s  windspeed  approximated  the  average  ascent  velocity  of  a  balloon 
which  typically  ranges  from  4  to  6  m/s.  While  it  was  in  the  wind  tunnel. 
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we  illuminated  the  probe  with  an  unfocused  visible  laser  diode  to  provide  a 
periodic  instantaneous  heating  effect.  This  was  to  simulate  the  effect  of  an 
instantaneous  change  in  air  temperature  as  the  probe  moved  through  different 
air  masses  as  it  would  throughout  a  flight.  The  control  portion  of  this  experiment 
was  conducted  at  sea  level  conditions.  The  wind  tunnel  was  built  of  styrofoam 
in  order  to  minimize  the  influence  of  any  external  heat  sources  on  the 
thermocouple.  Figure  24  Illustrates  the  experimental  setup.  Appendix  C  contains 
additional  details  of  the  experiment. 
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Figur*  24.  Schematic  of  wind  Tunnel  Thermocouple  Response 
Experiment  Conducted  in  Monterey,  CA  at  0  km  altitude  and  25* 
C  on  8  November  1992. 
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3.  Experimental  Results 

The  12.4  pm  probe  had  a  measured  time  constant  of  5.6  +/-  0.2  ms 
in  a  4.4  m/s  wind  in  the  wind  tunnel  at  0  km  altitude.  Figure  25  is  a  plot  of  the 
probe  response  to  the  periodic  laser  puise  at  this  windspeed.  It  illustrates  the 
rapid,  but  not  instantaneous,  response  of  the  probe  to  an  instantaneous 
temperature  change.  The  horizontal  distance  along  the  x-axis  between  the 
instantaneous  change  and  the  time  it  takes  the  amplitude  to  reach  the  1-1 /e  point 
is  the  time  constant  r.  Appendix  C  contains  plots  of  the  results  of  performing  this 
test  on  a  25  pm  copper-constantan  probe,  a  probe  used  frequently  in  previous 
experiments.  Table  X  summarizes  these  data. 


WIND  TUNNEL  PROBE  RESPONSE 

Typ*  "E"  1  2.4  micron  Prob*  (V=s4.4  m/s) 


Rgure  25.  Experimental  Response  of  a  12.4  pm  Diameter  Thermocouple  Probe 
Illuminated  by  a  Periodic  Laser  Diode  in  a  4.4  m/s  Wnd  in  the  NPS  Wind  Tunnel 
Experiment  of  2  November  1992.  Time  Constant  Measured  to  be  r  =  5.6  +/-  0.2 
ms. 
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TABLE  X.  Measured  Probe  Response  in  Wind  Tunnel  for  12.4  and  25  //m 
Wires. 


Wind 

(m/s) 

12.4  fjm 
Chrom-Const 

25  /jm  Copper- 
Const 

Response  Ratio 
(25;/m/12.4//m) 

0 

15.7 

69 

4.4 

1.5 

7.1 

28.5 

4.0 

3.0 

6.4 

21 

3.3 

4.4 

5.6 

19.5 

3.5 

These  results  show  a  factor  of  four  increase  in  probe  response  time 
for  a  factor  of  two  increase  in  probe  diameter.  This  shows  the  effect  of  the 
thermal  mass  slowing  the  response  time.  The  quadratic  diameter  dependence 
arises  from  Eq.  (21)  because  the  h^  term  depends  on  d.  Rearranging  Eqs.  (16), 
(17),  and  (20)  yields  an  equation  of  the  form  C  x  (dVd")  where  C  is  a  grouping 
of  constants  and  n  is  an  empirical  constant  equal  to  0.33  for  the  Reynolds 
number  regime  of  the  12.4  ^m  probe  with  a  velocity  of  4.4  m/s.  The  d®/d"  term 
then  reduces  to  d^  which  accounts  for  the  approximate  d“  ratio  between  the 
increase  in  the  time  constant  based  on  a  given  change  in  probe  diameter. 
Applying  this  same  rationale  to  the  spherical  junction  beads  for  the  two  different 
probe  sizes  should  account  for  most  of  tfie  remaining  difference. 

4.  Comparison  of  Experiment  vs.  Theoretical  Calculations 

As  previously  shown,  Is  a  complicated  function  of  numerous 
variables;  therefore  the  characteristic  time  constant,  r,  is  also  a  complicated 
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function  that  should  change  in  a  similar  fashion.  Table  XI  compares  experimental 
response  measurements  with  analytical  model  calculations  for  the  12.4  //m  probe. 
Appendix  C  contains  the  experimental  data  and  the  results  of  the  regression 
analysis.  We  did  not  include  the  plots  for  the  11.4  and  13.4  /ym  probes  here 
because  they  would  obscure  the  results  and  because  their  relationship  to  the 
response  of  the  12.4  /ym  probe  is  shown  in  Figures  22  and  23.  Figure  26  plots 
the  time  constant  at  0  km  as  a  function  of  windspeed.  This  shows  that  at  v  > 
1.5  m/s  windspeed,  the  analytical  model  closely  approaches  the  wind  tunnel 
experimental  measurements  for  a  12.4  //m  wire  with  a  26  /;m  bead.  At  v  =  4.4 
m/s,  the  average  ascent  velocity,  the  curves  of  the  experimental  regression  data 
virtually  overlay  the  modeled  calculations.  At  v  ^  1.5  m/s  the  plot  shows  the 
damping  effect  of  the  larger  thermal  mass.  Rgure  27  shows  the  same  data  using 
Eq.  (24),  the  anemometer  response  expression. 

TABLE  XI.  Time  Constant  Comparison;  Experimental  Measurements  vs. 
Analytical  Model  for  the  NFS  12.4  +/-  1.0//m  Chromel-Constantan  Thermocouple 
Probe. 


1 . 

V 

EXPERIMENT 

EXPERIMENT 

MODEL 

MODEL 

(m/s) 

Wind  Tunnel 

Curve  fit 

12.4  fjm 

12.4  ^/m 

No  bead 

26  ^m  bead 

0.01 

15.7 

15.0 

14.7 

19.4 

1.5 

7.1 

7.5 

6.0 

7.9 

3.0 

6.4 

6.2 

4.8 

6.3 

m 

5.6 

5.6 

mmm 

5.6 

*  Average  ascent  velocity 
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FIgurs  26.  Thermocouple  Response  Times  vs.  Windspeed.  Comparison  of 
Experimental  (Wind  Tunnel)  Measurements  with  Analytical  Model  Calculations  for 
12.4  +/- 1 .0  //m  Chromel-Constantan  Wire,  with  and  without  Correction  for  26  //m 


Junction  Bead. 
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Rgura  27.  1/e  Thermocouple  Probe  Response  vs.  Square  Root  of  Wind  Velocity. 
Comparison  of  Wind  Tunnel  Measurements  with  Analytical  Model  for  12.4  +/- 1.0 
fim  Chromel-Constantan  Wire,  with  and  without  Correction  for  26  um  Junction 


Bead. 
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For  a  12.4  //m  probe,  the  analytical  model  overestimates  the  observed 
probe  response  by  29%  when  there  is  no  correction  factor  for  the  bead.  When 
a  correction  for  the  bead,  assumed  to  be  spherical  and  with  a  diameter  of  26  //m, 
is  made,  the  model  is  within  1%  of  the  experimental  probe  response.  It  is  not 
probable  that  the  one  correction  made  for  bead  diameter  was  sufficient  to 
improve  the  model’s  accuracy  by  a  factor  of  30.  Wind  tunnel  experimental  biases 
such  as  uneven  heating  conditions,  non-uniform  windspeed,  or  other  slight 
deviations  from  standard  atmospheric  conditions  in  the  laboratory  during 
measurements  could  have  resulted  in  a  larger  difference  between  the  model  and 
experimental  results.  Coupled  with  this  is  the  uncertainty  in  the  exact  dimensions 
of  both  the  thermocouple  wire  and  junction  bead  diameters.  This  uncertainty, 
although  not  insignificant,  does  not  obscure  the  actual  decrease  in 
responsiveness,  by  a  factor  of  approwmately  4  (corresponding  to  the  decrease 
in  hg  by  the  same  factor),  that  occurs  as  the  probe  ascends  from  0  to  30  km  as 
shown  in  Figures  25  and  26.  Figures  22  and  23  show  that  variations  in  the  wire 
diameter,  even  within  the  manufacturer's  specified  tolerances,  can  alter  response 
times  by  22%  under  experimental  conditions  which  equates  to  2  ms  at  10  km  and 
3  ms  at  20  km.  Table  XII  summarizes  these  data. 
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TABUE  Xli.  Differences  in  Analytical 
Model  Probe  Response  due  to 
Tolerance  in  Wire  Diameters  as  a 
Function  of  Altitude  for  the  NPS  12.4  +/- 
1 .0  fjm  Thermocouple  Probe. 


ALT 

(km) 

Min  r 
(ms) 

Max  r 
(ms) 

Lt  I 
(ms)  j 

0 

n 

6.2 

1.2  1 

10 

8.4 

10.4 

1 

14.3 

IT’S 

3.3  1 

1  30 

23.3 

28.8 

5.5  1 

D.  PROBE  RESPONSE  SUMMARY 

The  tolerance  in  the  diameter  of  the  12.4  //m  Type  “E"  chromel-constantan 
thermocouple  wire  and  junction  bead  make  analytical  predictions  of  its  true 
response  uncertain  to  ^proximately  11%  (5.6  +/-  0.6  ms  at  0  km).  Additionally, 
the  3:1  ratio  of  the  bead  to  wire  diameter  results  in  higher-order  response  effects. 
The  12.4  >t/m  chromel-constantan  probe  used  by  NPS  is  four  times  faster  than  the 
25//m  copper-consteuitan  probe  used  extensively  in  previous  experiments  due  to 
its  relative  mass  to  surface  area  ratio. 

The  effect  of  the  spherical  weld  junction  on  the  thermocouple  wire  was  to 
slow  the  response  time  by  30%.  Correcting  for  a  spherical  junction  bead  allowed 
the  12.4  //m  probe  to  achieve  a  modeled  performance  within  1 1  %  of  experimental 
wind  tunnel  results.  A  probe  with  uniform  diameter  and  a  sphericed  bead,  or  no 
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bead,  would  allow  for  even  more  precise  extraqpolations  of  modeled  response 
calculations.  A  faster  probe  response  would  improve  the  vertical  layer  resolution 
as  well  as  permit  more  precise  structure  function  measurements  at  smaller  probe 
spacings.  A  probe  that  optimizes  response  by  minimizing  its  wire  diameter,  weld 
bead  diameter,  density,  and/or  its  specific  heat  capacity  would  achieve  this  result. 
A  thermocouple  made  of  the  10  //m  wire  currently  manufactured  by 
Isabellenheutte  Heusler,  Gmbh,  in  Germany  should  be  60%  faster  than  the  12.4 
jwm  chromel-constanteui  probes  used  by  NPS,  assuming  a  proportionately  sized 
junction  bead  of  ~21  //m. 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

This  thesis  investigated  two  systematic  error  effects  in  balloon  theimosonde 
systems  used  by  both  NPS  and  the  U.S.  Air  Force  Geophysics  Directorate,  it 
supports  the  existence  of  an  order  of  magnitude  solar  contamination  of  all 
thermal  probe  daytime  data  above  the  tropopause.  This  most  probably  arises 
from  a  combination  of  package  motion  (rotation  or  pendulum-like  oscillations 
producing  a  low  frequency  signal),  and  differential  solar  heating  of  the  probes 
due  to  probe  wire  osdilation.  The  high  correlation,  91%,  observed  for  the 
temperature  differences  in  the  nighttime  flight  and  the  low  correlation  obtained 
during  ail  daytime  flights,  ranging  from  50%  to  62%,  particularly  near  the 
tropopause,  support  the  hypothesis  that  the  apparent  order  of  magnitude 
increase  in  optical  turbulence  is  a  solar  induced  instrument  effect. 

This  thesis  also  investigated  the  tiiermal  response  time  of  the  thermocouple 
used  in  the  NPS  differential  thermosonde.  The  relatively  large  size  of  the 
thermocouple  made  it  imperative  to  reconcile  its  modeled  response  as  a  function 
of  altitude  witii  controlled  laboratory  tests.  The  inconsistencies  in  the  diameter 
of  the  12.4  ^m  type  "E*  chromel-constantan  thermocouple  wire  and  its  spherical 
junction  bead  limit  predictions  of  its  true  response  to  a  precision  of  11%. 
Additionally,  the  6:1  ratio  of  the  bead  to  wire  diameter  results  in  higher-order 
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response  effects.  Based  upon  calculated  changes  in  atmospheric  heat 
conductivity  the  thermal  response  time  for  the  12.4  //m  chromel-constantan  probe 
should  be  approximately  9-10  ms  at  sdtitudes  of  10  km  and  roughly  15-17  ms  at 
20  km. 

The  geometric  characteristics  of  the  wire  and  the  junction  bead  that  govern 
probe  response  can  vary  for  each  probe.  Measuring  the  time  constant  for  a 
probe  in  a  wind  tunnel  experiment,  and  then  using  the  same  probe  in  a  launch, 
would  have  the  effect  of  calibrating  that  probe  against  the  model. 

B.  RECOMMENDATIONS 

To  eliminate  solar  effects  apparent  in  the  daytime  launches  due  to  rotation 
and  oscillation  of  the  sonde,  it  is  necessary  to  stabilize  the  sonde  as  much  as 
possible.  The  most  reliable  and  consistent  data  will  occur  when  the  probes 
maintain  a  consistent  alignment  with  respect  to  the  sun.  Possible  solutions  to 
this  may  include  the  use  of  a  gyro  mechanism  to  stabilize  the  platform,  vertical 
baffles  or  fin  stabilizers,  solar  celts  to  determine  relative  solar  angle,  and  weights 
or  anchors  to  hold  the  apparatus  steady.  To  reduce  the  effects  of  probe  wire 
vibration,  the  tungsten  wire  probes  must  be  attached  straight  across  the  probe 
supports  rather  than  in  a  flexible  configuration.  This  will  involve  a  delicate 
compromise  to  avoid  thermal  stress  and  strain  on  the  wire  and  its  contacts  and 
yet  suppress  any  vibration. 
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Comparison  of  results  from  taut  wire  probes  with  those  from  the  GP-type 
loose  wire  thermocouple  should  help  clarify  the  issue  of  probe  vibration. 
Calibrating  a  set  of  probes  in  a  wind  tunnel  experiment  prior  to  a  launch  should 
improve  consistency  in  thermosonde  temperature  measurements, 

A  faster  probe  would  improve  measurement  precision  of  the  NPS  system 
in  frie  stratosphere.  A  probe  that  optimizes  response  by  minimizing  its  size, 
density,  and/or  its  specific  heat  capacity  would  achieve  this  result.  Using  a  10 
fjm  thermocouple  probe  wire,  available  in  Germany,  should  reduce  the  response 
time  by  a  factor  of  about  1 .6.  Additionally,  ensuring  that  thermocouple  wire 
junction  beads  are  geometrically  uniform  would  allow  for  more  precise 
extrapolation  of  experimental  data  to  stratospheric  data  Regeu’diess  of  the 
probes  used,  thorough  consistency  checks  should  be  built  into  subsequent 
experiments.  Thermocouple  probes  from  the  same  lot  should  be  measured  for 
wire  diameter  and  bofri  bead  size  and  shape  under  a  calibrated  SEM.  These 
same  probes  should  then  be  tested  for  response  in  a  wind  tunnel  experiment, 
varying  the  laser  to  illuminate  botti  the  wires  and  the  junction  beads.  Finally, 
these  probes  should  be  launched  In  a  thermosonde  package  that  has  redundant 
means  of  measurement  at  equal  distaices  to  establish  an  in-flight  empirical 
database  with  which  to  resolve  atmospheric  structure  functions. 
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APPENDIX  A 


DIURNAL  VARIATION  EXPERIMENT  2  DATA 

Additional  data  from  OVE  2  are  shown  below.  The  data  contained  in  this 
appendix  were  used  in  the  statistical  summary  shown  in  Table  11  and  do  not 
exhibit  any  atypiced  behavior.  Rgures  28  thru  33  show  the  results  of  the  4,  5,  and 
6  December  1990  flights.  Rgures  34  thru  37  show  the  averaged  profiles  for 
the  7  and  8  December  flights. 
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Figure  28.  Log  (Raw  C„^Profiie  (Day),  Collected  by  the  Geophysics  Directorate, 
from  1.2  to  30  km,  Holloman  Air  Force  Base,  NM,  4  December  1990  (Probe  Set 
1). 
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Rgure  29.  Log  (Raw  C„*)  Profile  (Day),  Collected  by  the  Geophysics  Directorate, 
from  1,2  to  30  km,  Holloman  Air  Force  Base,  NM,  4  December  1990  (Probe  Set 
2). 


Ffgur*  30.  Log  (Raw  C„*)  Profile  (Day),  Collected  by  the  Geophysics  Directorate, 
from  1.2  to  30  km,  Holloman  Air  Force  Base,  NM,  5  December  1990  (Probe  Set 
1). 


Fi9ur«31.  Log  (Raw  Profile  (Day),  Collected  by  the  Geophysics  Directorate, 
from  1,2  to  30  km,  Holloman  Air  Force  Base,  NM,  5  December  1990  (Probe  Set 
2). 


Rgur*  32.  Log  (Raw  C„^  Profile  (Day),  Collected  by  the  Geophysics  Directorate, 
from  1,2  to  30  km,  Holloman  Air  Force  Base,  NM,  6  December  1990  (Probe  Set 

11. 
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Figure  33.  Log  (Raw  C„*)  Profile  (Day),  Collected  by  the  Geophysics  Directorate, 
from  1.2  to  30  km,  Holloman  Air  Force  Base,  NM,  6  December  1990  (Probe  Set 
2). 


85 


86 


Figure  35.  Averaged  Profile  (Day),  Collected  by  the  Geophysics  Directorate, 
from  1.2  to  30  km,  Holloman  Air  Force  Base,  NM,  7  December  1990  (Probe  Set 
2). 
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Rgur«  36.  Averaged  Profile  (Night),  Collected  by  the  Geophysics 
Directorate,  from  1.2  to  30  km,  Holloman  Air  Force  Base,  NM,  8  December  1990 
(Probe  Set  1). 


Rgur«  37.  Averaged  C„*  Profile  (Night),  Collected  by  the  Geophysics 
Directorate,  from  1.2  to  30  km,  Holloman  Air  Force  Base,  NM,  8  December  1990 
(Probe  Set  2). 


APPENDIX  B 


THERMOCOUPLE  PROPERTIES 

Table  XIII  shows  the  specific  heat  capacities  of  the  materials  that 
comprise  the  type  E  chromel-constantan  probes  [Ref.  18:p.  12-107],  Chromel- 
Constantan  is  75%  chrome!  and  25%  constantan.  Chromel  is  90%  nickel  and 
10%  chromium.  Constantan  is  57%  copper  and  43%  nickel  [Ref.  27:p.  133], 
Applying  these  percentages  to  the  specific  heat  capacities  in  Table  XIII  results 
in  a  specific  heat  capacity  of  C =0.4367  (J/g*K)  for  the  chromel-constantan 
thermocouple  probe. 


TABLE  XIII.  Type  E  Chromel-Constantan  Thermocouple 
Material  Specific  Heat  Capacity,  C. 


Nickel 

(Ni) 

C =0.1 06 

Chromium 

(Cr) 

C=0.110 

Copper 

(Cu) 

C=0.092 

Chromel 

90% 

10% 

Constantan 

43% 

57% 

Specific  heat  capacity,  C,  is  in  (J/g*K). 
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APPENDIX  C 


WIND  TUNNEL  PROBE  RESPONSE  EXPERIMENT 


Table  XIV  contains  the  wind  tunnel  measurements  on  the  12.4  //m  probe. 

TABLE  XIV.  Wind  Tunnel  Probe  Response  Measurements  vs.  \Afindspeed, 
12.4  jjm  Chromel-Constantan  Probe,  8  November  1992. 


f  V  (m/s) 

Measured  r 
(ms) 

SORT  of  V 

1/r 

0 

15.7 

0 

0.06 

1.5 

7.1 

1.22 

0.14 

3.0 

6.4 

1.73 

0.16 

4.4 

5.6 

2.1 

0.18 

Regression  analysis  on  this  data  yielded  values  of  a  -  0.066  and  b  0.054 
for  the  constants  used  in  Eq.  (24). 


Figure  38  plots  the  response  of  a  25  ;/m  copper-constantan  probe  in  the 
wind  tunnel  in  which  the  time  constant  was  measured  to  be  r  =  19.5  +/-  1.0 
ms. 


25  Micron  Chromel-Constantan  Probe 

3  m/s  Wfndspeed 


Figure  38.  Experimental  Thermocouple  Response  of  a  25  //m  Copper- 
Constantan  Probe  illuminated  by  a  Periodic  Laser  Diode  Measured  in  a  Wind 
Tunnel  at  v  =  4.4  m/s  in  the  NPS  Wind  Tunnel  Experiment  of  8  November  1992. 
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Table  XV  contains  the  instrument  settings  for  each  item  of  equipment 


used  in  the  wind  tunnel  probe  response  experiment. 


TABLE  XV.  Wind  Tunnel  Experiment  Instrument  Settings. 


NICOLET  PRO  30 
DIGITAL 

OSCILLOSCOPE 

HP3314A 

FUNCTION 

GENERATOR 

HP6236B 

POWER 

SUPPLY 

LASER 

DIODE 

2.02  hz 

0-24  V 

-  15  cm 
from  probe 

500  sweeps 

4.26  V  Ampl. 

120  mv 

2.02  V  Offset 

_ J 

93 


LIST  OF  REFERENCES 


1.  Hufnagel,  R,E.,  "Propagation  Through  Atmospheric  Turbulence,"  in  The 
Infrared  Handbook.  The  Infrared  Information  Analysis  Center,  Environmental 
Research  Institute  of  Michigan,  19^. 

2.  Brown,  J.H.  and  others,  Study  of  Possible  Solar  Heating  Effects-Error 
Analysis.  GL-TR-89-0178,  Geophysics  Laboratory,  Hanscom  Air  Force 
Base,  Massachusetts,  July  1989. 

3.  Weitekamp,  M.  R.,  Investigation  of  the  Source  of  Thermosonde  Measured 
DiumalVariation  of  Optical  Turbulence.  M.S.  Thesis,  NavaJ  Postgraduate 
School,  Monterey,  California,  December  1990. 

4.  Gast,  V.J.,  Characterization  of  Atmospheric  Turbulence  For  High  Resolution 
Imaging  and  Laser  Propagation  ,  M.S.  Thesis,  Naval  Postgraduate 
School,  Monterey,  California,  December  1992. 

5.  Wilson,  J.  and  Hawkes,  J.F.B.,  Optoelectronics.  An  Introduction. 

Prentice  Hall,  1989. 

6.  Krieth,  Principles  of  Heat  Transfer.  International  Textbook  Company,  1968. 

7.  White,  F.M.,  Fluid  Mechanics.  New  York,  McGraw  Hill,  1979. 

8.  Tatarski,  V.l.,  Wave  Propagation  in  a  Turbulent  Medium.  McGraw  Hill,  1961. 

9.  Walters,  D.  L.  and  Kunkel,  K.  E.,  "Atmospheric  Modulation  Transfer 
Function  For  Desert  and  Mountain  Locations:  the  Atmospheric  Effects  on  r^". 
Journal  of  the  Optical  Society  of  America.  Volume  71,  No.  4,  April  1981. 

10.  Ochs,  G.R.  and  others,  "Laser-Beam  Scintillation  over  Horizontal  Paths  from 
5.5  to  145  Wlometers,"  In  Journal  of  the  Optical  Society  of  America.  Vol.  59, 
No.  2,  February  1969. 


94 


11.  Bauer,  E,  Atmospheric  Effects  on  Airborne  Lasers  (ABU  for  Tactical 
Missile  Defense  (TMD1:  Clouds  and  Turbulence  Draft  Report.  IDA  Log 
No.  HQ  91-40634/1,  Institute  for  Defense  Analyses,  Alexandria,  Virginia, 
December  1991. 

12.  Brown,  J.H.  and  others,  Sonde  Experiments  for  Comparative 
Measurements  of  Optical  Turbulence.  AFGL-TR-82-0079,  ADA  1 1 8740,  24 
February  1982. 

13.  Hummel,  J.R.  and  Shettle,  E.P.,  Effects  of  Solar  Heating  bv  Aerosols  and 
Trace  Gases  on  the  Temperature  Structure  Constant.  Scientific  Report 
LTR90-009  written  by  Sparta,  Inc.  for  Geophysics  Laboratory,  Hanscom  AFB, 
Massachusetts  (Report  GL-TR-90-0349),  9  August  1990. 

14.  Beland,  R.,  Experiment  Plan  for  DVE2  f Diurnal  Variation  Experiment  2). 
internal  Geophysics  Laboratory  memorandum,  23  November  1990. 

15.  Kiefer,  D.A.,  Development  of  a  Single  Point  Temperature  Probe  to 
Measure  the  Temperature  Structure  Parameter.  M.S.  Thesis,  Naval 
Postgraduate  School,  Monterey,  California,  June  1990. 

16.  Holman,  J.P.,  Heat  Transfer.  McGraw  Hill  Book  Company,  New  York, 

1990. 

17.  CRC  Handbook  of  Chemistry  and  Physics,  ed.  Lide,  David  R.,  71st  Ed., 

1991. 

18.  OMEGA  Complete  Temperature  Measurement  Handbook  and 
Encyclopedia,  vol.26.  Omega  Engineering,  Stamford,  Connecticut,  1988. 

19.  Handbook  of  Geophysics  and  the  Space  Environment,  ed.  A.S.  Jursa,  Air 
Force  Geophysics  Laboratory,  19K. 

20.  Telephone  conversation  between  John  S.  Daley,  of  Omega  Engineering,  Inc. , 
Quotations  Engineering  Department,  and  the  author,  3  December  1992. 

21.  Kothandaraman,  C.P.  and  Subramanyan,  S.,  Heat  and  Mass  Transfer  Data 
Book.  Halsted  Press,  New  York,  1975. 

22.  Lunardini,  V.J.,  Heat  Transfer  in  Cold  Climates.  Van  Nostrand  Reinhold 
Company,  New  York,  1981. 


95 


23.  Kacak,  S.  and  others,  Handbook  of  Single-Phase  Convective  Heat 
Transfer.  John  VWley  and  Sons,  New  York,  1987. 


24.  Moffat,  R.J.,  "Gas  Temperature  Measurement,"  in  Temoerature-lts 
Measurement  and  Control  in  Science  and  Industry.  Vol.  3,  Part  2,  Rheinhold 
Publishing  Company,  New  York,  1962. 

25.  StereoScan  200,  TL2025-OM,  Issue  3,  April  85,  Cambridge,  England. 

26.  Hoover,  C.R.,  Investigations  of  a  Single  Point  Temperature  Probe  for 
Measurement  of  Atmospheric  Turbulence.  M.S.  Thesis,  Naval  Postgraduate 
School,  Monterey,  California,  December  1991. 

27.  Wnzie,  P.A.,  Thermocouple  Temperature  Measurement.  New  York,  John 
Wiley  and  Sons,  1973. 


96 


INITIAL  DISTRIBUTION  UST 


1.  Defense  Technical  Information  Center  2 

Cameron  Station 

Alexandria,  VA  22304-6145 

2.  Library,  Code  52  2 

Naval  Postgraduate  School 

Monterey,  CA  93943-5002 

3.  Professor  Karlheinz  E.  Woehler  1 

Chairman,  Department  of  Physics 

Naval  Postgraduate  School  (Code  PHWh) 

Monterey,  CA  93943 

4.  Professor  Donald  L  Walters  6 

Naval  Postgraduate  School  (Code  PHWe) 

Monterey,  CA  93943 

5.  U.S.  Air  Force  Geophysics  Directorate  2 

Phillips  Lab  (GPOA) 

ATTN;  Bob  Beland 
Hanscom  AFB,  MA  01731-5000 

6.  DCSOPS  2 

Space  and  Special  Weapons  Directorate 

ATTN:  DAMO-SWS  (CPT  Roper) 

Room  2C534,  The  Pentagon 
Washington,  DC  20310 


97 


